
Figure 1. 



Table 1. Properties of galaxies in the catalog. 
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-37 
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1344 


287 
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589 
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135 
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67.1 
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11.48 
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56.7 
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3344 
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424 
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3351 
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777 
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47.5 
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46.2 
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-5 


10:47:09.1 


+ 14:02:35 
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691 
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3379 
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10.24 


888 
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3486 


Lowell 
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424 


42.2 
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Lowell 


6 
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693 
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11:12:27.9 
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11.95 


1191 
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+ 13:09:31 


10.25 


806 
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72.8 
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3631 


Lowell 
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11:18:13.3 
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11.01 


1157 
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17.3 
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-09:44:33 
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62.1 


12 
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765 
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58.3 
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+47:05:49 
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46.2 
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Table 1. (Continued.) 



NGC 


Obs. 


T 


a 


(1950) 


h (1950) 




{v)hc 


D25 


/ 


p. a. 


41zd 


T n 

Loweii 


5 


12: 


:07:30.9 


1 r\o r ^ 1 ^ 

+02:57:17 


11, 


.98 


1 ooo 

1328 


O/^ 1 

261 


/I 

42.2 


1 r 

135 


4125 


T n 

Loweii 


-5 


12: 


:07:25.1 


1 n ^ 1 o /II 

+65:13:41 


10, 


.65 


1 o r /I 

1354 


/I r 

345 


56.7 


f\ r 

95 


4136 


T n 

Loweii 


5 


12: 


:06:45.7 


1 o r\ 1 o 1 o 

+30:12:18 


11, 


.69 


n r\'~7 

607 


000 

238 


01 1 

21.1 





4144 


T n 

Loweii 


6 


12: 


:07:28.3 


1 An A A r\'i 

+46:44:07 


12, 


.05 


f^n '~7 

267 


0/^1 

361 


'I'l A 

77 A 


104 


4157 


T n 

Loweii 


5 


12: 


:10:26.7 


1 o o or\ 

+50:33:20 


12, 


.66 


oor\ A 

2294 


1 /I 

134 


r\ 1 

39.1 


1 1 r 

115 


A 1 

4178 


Jr^aiomar 


8 


12: 


:10:13.8 


1 11 r\Oi A Oi 

+ 11:08:48 


11, 


.90 


o ''vn 

376 


307 


n f\ 

69.2 


r\ 

30 


4189 


Jr^aiomar 


6 


12: 


:11:13.9 


1 1 O /I O 11 

+ 13:42:11 


12, 


.51 


Olio 

2112 


1 /I 

143 


/I 

43.6 


r 

85 


4192 


Jr^aiomar 


2 


12: 


11 1 r\ 

:11:15.9 


1 1 1 r\ A f\ 

+ 15:10:49 


10, 


.95 


1 /I o 

-142 


n 

586 


^0 

73.6 


1 r r 

155 


4216 


Jr^aiomar 


3 


12: 


:13:20.9 


1 1 o or OO 

+ 13:25:22 


10, 


.99 


1 or\ 

129 


487 


77 A 


1 f\ 

19 


4242 


T n 

Loweii 


o 

8 


12: 


:16:57.5 


1 AT" A C\ ^'1 

+45:40:37 


11, 


.37 


r 1 ^ 

517 


300 


A C\ '1 

40.7 


r 

25 


4254 


Jr^aiomar 


5 


12: 


1 1 o r\ 

:16:18.0 


\ -\ A A -\ A C\ 

+ 14:41:42 


10, 


.44 


o /I r\ ^ 

2407 


000 

322 


or\ A 

29.4 





4258 


Jr^aiomar 


4 


12: 
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-5 
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4394 


Jr^aiomar 
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217 
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Loweii 
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12: 
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314 


/I 
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.36 


OO 

2278 


/I 1 r 
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4526 
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-2 
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:33:34.0 
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70.7 
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66.0 
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4564 
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-5 


12: 


:35:51.2 
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12, 
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65.4 


47 


4569 
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2 
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:34:17.9 
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23 
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Palomar 
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Table 1. (Continued.) 
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80 
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Jr^aiomar 


6 
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293 


76.1 
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35 
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6 
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.90 
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95 
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Loweii 


3 


12:51:49.8 
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11, 
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57.5 
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23 


4826 
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Loweii 


2 


12:56:08.7 
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9, 


.36 


/I 1 o 

412 
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r ^ r 

57.5 
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115 


4861 
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Loweii 


9 


12:56:40.3 


+35:07:56 


12, 


.90 
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843 
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238 
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15 


4866 


T n 
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-1 


12:56:57.9 
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12, 
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O 
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77.7 


87 
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Loweii 


4 
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+37:06:53 
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65 
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353 
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48.6 


f\ r 

95 
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Loweii 


5 


13:50:20.3 


-00:52:05 


11, 


.99 


1 o o o 
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43.6 
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15 
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50 


5701 


T n 

Loweii 





14:36:41.5 
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.45 
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145 
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5850 
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Loweii 


3 
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15:06:39.4 
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13 
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Lowell 


6 
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+62:20:17 


11, 
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66.5 


28 
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Lowell 


6 
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Table 2. FITS header keys. 



N 



ame 



Type Description 



EXPOSURE float Length of exposure in seconds. 

FILTERl character Photometric passband used for the observation. 

RA character Right ascension of the pointing center during the observation. 

DEC character Dechnation of the pointing center during the observation. 

AIRMASS float Airmass at the beginning of the exposure. 

TIME character Universal time (UT) at the begining of the exposure. 

SKY float The background sky level (in ADU). 

SKYSIG float RMS noise per pixel of the sky (in ADU). 

DNATO_ST float Counts (in ADU) representing magnitude in the image, 

calibrated using observations of standard stars. 

DNATO_BV float Counts (in ADU) representing magnitude in the image, 

calibrated using B and B — V from RC3. 

G_CENT_X integer Horizontal (x) coordinate of the photometric center of the 

galaxy in the trimmed image, in pixels. 

G_CENT_Y integer Vertical (y) coordinate of the photometric center of the 

galaxy in the trimmed image, in pixels. 

SATURAT integer Equals 1 if the galaxy center was not saturated, 2 if it was 

saturated (and divided by 2). 

PSF_FWHM float Width (in arcsec) of the point spread function, obtained 

from stars before foreground star removal. 

B_RC3 float Total (asymptotic) Johnson B magnitude from RC3, 

column 6, line 1. 

B-V_RC3 float Total (asymptotic) Johnson B — V color from RC3, 

column 7, line 1. 

PA_RC3 float Major axis position angle, (measured from North eastward, 

between 0° and 180°) from RC3, column 5, line 1. 

B0A_RC3 float Decimal logarithm of the ratio of the major axis isophotal 

diameter to the minor isophotal diameter, both reduced to 
jiB = 25 mag arcsec"^, from RC3, column 4, line 2. 

SIZE_RC3 float Decimal logarithm of the apparent major axis isophotal 

diameter reduced to jis = 25 mag arcsec"^, 
from RC3, column 4, line 1. 

TYPE_RC3 character Numerical morphological type T from RC3, column 3, line 1. 

VEL0_RC3 float Weighted mean heliocentric radial velocity, calculated 

from RC3, column 10, line 4. 




By. (mag) 



Figure 3. 




Figure 4. 




Figure 5. 



A CATALOG OF DIGITAL IMAGES 
OF 113 NEARBY GALAXIES 



ZSOLT FREI, 
PURAGRA GUHATHAKURTA^ 
AND 

James E. Gunn 

Department of Astrophysical Sciences 
Princeton University, Peyton Hall 
Princeton, NJ 08544 

AND 

J. Anthony Tyson 

AT&T Bell Laboratories 
Murray Hill, NJ 07974 

Submitted to the Astronomical Journal 

Received 1995 Accepted 199 

^ Current address: UCO/Lick Observatory, University of California, Santa C 
CA 95064 



-2- 



ABSTRACT 



We present a digital catalog of images of 113 galaxies in this paper. These 
galaxies are all nearby, bright, large and well resolved. All images were recorded 
with charge coupled devices (CCDs) at the Palomar Observatory with the 1.5 meter 
telescope and at the Lowell Observatory with the 1.1 meter telescope. At Palomar 
we used the Thuan-Gunn </, r and i photometric bands (Thuan and Gunn, 1976) to 
take 3 images each of 31 spiral galaxies; at Lowell we used the B j and R bands (2 
images per galaxy) of the photometric system by GuUixson eA al. (1994) to observe 
82 spirals and ellipticals. The galaxies were selected to span the Hubble classification 
classes. All data are photometrically calibrated with foreground stars removed. 
Important data on these galaxies published in the "Third Reference Catalog of 
Bright Galaxies" (de Vaucouleurs eA al., 1991) are recorded in the FITS file headers. 
All files are available through anonymous FTP from "astro.princeton.edu", and 
Princeton University Press will soon publish the data on CD-ROM. 
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1. INTRODUCTION 

We describe a digital catalog of galaxy images. The direct motivation for 
the development of such a catalog was the need for a suitable data set to test the 
automatic galaxy classification techniques we are currently working on (Frei, 1995a, 
in preparation). Automatic means of classifying galaxies are necessary to handle 
the huge amounts of imaging data that we soon to be available from large survey 
projects, such as the Sloan Digital Sky Survey (Gunn and Knapp, 1993). We think, 
however, that these images are useful beyond the scope of this project; consequently, 
we make the digital images available. The catalog contains FITS files (Wells et al., 
1981) of the images, with an extensive set of corresponding information in the file 
headers. 

To test the classification techniques, we needed a set of galaxies that span 
the full range of the Hubble morphological classes. Another requirement was to 
have large, well resolved galaxies so that we could develop our techniques step-by- 
step: first, solve the easiest challenges that large galaxies pose, then improve the 
techniques so that they work for galaxies which are further away. The performance 
and limitations of the classification system can be tested on degraded images. If 
images of galaxies in the digital catalog are available in several passbands then one 
can interpolate or extrapolate in color, thus synthetically "redshifting" the image of 
a given galaxy to any desired z (while decreasing the resolution at the same time). 

Consequently, we developed a catalog of images of nearby, well resolved galaxies 
that were recorded in several passbands with charge coupled devices (CCDs). 
Although CCDs have been in extensive use for about 15 years, there is no similar 
data set publicly available which we know of. 
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Early photometric surveys of galaxies date back to the 70's. Freeman (1970) 
collected surface photometric data for 36 spiral and lenticular galaxies to examine 
the disks. Kormendy (1977a, 1977b, 1977c) used photoelectric photometry to 
study a total of 20 red and neutral compact galaxies. King (1978) performed 
photographic surface photometry on 17 galaxies, mostly giant ellipticals, to compare 
the brightness profiles of giant ellipticals to those of star clusters in his famous study 
(King, 1962, 1965, 1966). Burstein (1979a, 1979b) investigated 18 SO galaxies. He 
also used photographic surface photometry to explore the luminosity profiles and the 
origin of SO galaxies. Boroson (1981) derived luminosity profiles from photographic 
photometry of 26 spiral galaxies. 

Kent (1984, 1985) did CCD surface photometry of field galaxies. His set 
contains 105 galaxies of different, widely distributed morphological types. These are 
not the largest galaxies on the sky; in fact, we found only three common between 
galaxies in his set and our set of 113 nearby galaxies. His aim in observing field 
galaxies was, among other things, to compare statistics of field galaxies to statistics 
of disk galaxies in clusters. Using the photometric data he calculated bulge-to-disk 
ratios based on intensity profiles along the major and minor axes. 

Other recent works on nearby galaxies with CCDs include Davis et al. (1985), 
who studied three well known elliptical galaxies in detail. Lauer (1985) had a 
larger sample of 42 galaxies, ellipticals and SOs, which allowed him to do systematic 
studies of central structures and core properties. Simien and de Vaucouleurs 
(1986) compiled a large data set to investigate bulge-to-disk ratios of spirals 
and lenticulars. Their data included 98 galaxies, 32 of which were observed by 
them. Other data came from works mentioned above (Boroson, 1981; Burstein, 
1979b; Kormendy, 1977b) and from two others (Yoshizawa and Wakamatsu, 1975; 
Whitmore and Kirshner, 1981); not all of these images were taken with CCDs. 
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Watanabe et al. (1982), Watanabe (1983), Okamura et al. (1984) and 
Watanabe eA al. (1985) studied more than 200 galaxies to develop a quantitative 
classification system. They used plates and digitized them with an isophotometer 
which used a one dimensional CCD as a detector. The culmination of their effort 
is the "Photometric Atlas of Northern Bright Galaxies" (Kodaira eA al., 1990). 
This book contains digitized images of 791 galaxies. Unfortunately, these data are 
available in a printed form only, and since the images are digitized from photographic 
plates, all problems presented by those plates (nonlinearity, for example) are 
incorporated into the digital form of the data. More recently, Han (1992) presented 
CCD surface photometry of spirals (284 galaxies in 16 clusters) and CoUess eA al. 
(1993) studied ellipticals and lenticulars (352 galaxies in distant clusters). However, 
the galaxies in the last two studies are too distant for the purposes of developing 
and testing our automated classification system. 

There are other catalogs with images of galaxies, all of which are printed 
from photographic plates (Sandage, 1961; Sandage and Bedke, 1988). Also worth 
mentioning are those efforts that resulted in a large collection of useful photometric 
data on galaxies, but which do not contain images (Sandage and Tammann, 1981; 
TuUy and Fisher, 1987; de Vaucouleurs eA al., 1991). 

Images of the galaxies in our catalog were obtained at two different sites. At 
the Lowell Observatory, we used the 1.1 meter telescope and a CCD camera. The 
photometric system is described in GuUixson eA al. (1994). We used two filters, 
B J and R. About 25% of our images were observed at the Palomar Observatory 
with the 1.5 meter telescope equipped with a CCD camera and the Thuan-Gunn 
photometric system (Thuan and Gunn, 1976). We used three passbands: g, r, and 
i. 

An unusual feature of our catalog is that foreground stars have been removed 
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from the frames. This has been done for several reasons: to derive quantitative 
parameters that describe the galaxy for the purposes of classification, we had to 
ensure that the foreground objects did not contribute to these parameters. Secondly, 
to synthetically "redshift" the images of these galaxies, the images had to be free of 
foreground stars. We think that star-removal will also be useful for other possible 
applications of this catalog. 

In the next section, we describe the observations. Section 3 contains details of 
the processing of the images. In Section 4, we tabulate and plot distributions of 
parameters of the galaxies in the catalog, and in Section 5 we describe the image file 
format and the information contained in the file headers. In the short discussion in 
Section 6, we give details on how to access the data through the computer network. 
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2. OBSERVATIONS 

2.1 Palomar 

Thirty-one galaxies were observed with the 1.5 meter telescope (P60) of 
Palomar Observatory on the night of 1991 May 4 with the Wide Field Prime 
Focus Universal Extragalactic Instrument (Wide Field PFUEI). Gunn and Westphal 
(1981) describe the original PFUEI, designed for the 5 meter telescope (P200) at 
Palomar. The wide field camera has a 306 mm f/4 aerial camera lens collimator and 
a Nikon 58 mm f/1.2 camera lens. The field of view is 16' X 16', and it is projected 
on an 800 X 800 Texas Instruments (TI) CCD in the Cassegrain focus, yielding a 
scale of l'!l9 pixel""^. The focal ratio of the system is f/1.65 (for more details see 
Hester and Kulkarni, 1989). 

Images were obtained in the </, r and i bands of the Thuan-Gunn photometric 
system. The original, 4-band system (u, g and r) was introduced by Thuan and 
Gunn (1976). It was later extended by Wade et al. (1979) with the i band, and by 
Schneider eA al. (1983) with the z band. Most CCDs are not efficient in the z band, 
but (/, r and i CCD surface photometry of galaxies have been extensively done in 
the past decade. These filters are centered at 500, 650 and 820 nm, respectively. 
See Frei and Gunn (1994) for a detailed description of the photometric system we 
used, including the system response curves (including the effects of atmospheric 
absorption, two reflections from aluminum, the transmission functions of fllters and 
optics, and the quantum efficiency of the detector). The 7.5 cm square fllters were 
placed in front of the camera optics to avoid problems with reflection. 

All exposures were 60 s. The 9 electron CCD readout noise was low enough 
to allow such a short exposure time and still ensure that the background was sky 



8 



noise limited. The gain (G) of the instrument is about 2 electrons/ ADU (analog- 
to-digital units). The effective noise in electrons (cTg) is the quadratic sum of the 
readout noise and the sky noise (the latter, ideally, being the square root of the 



where R is the readout noise in electrons and N is the counts in ADU in a given 
pixel. With 60 s exposure time, the level of the sky was typically around 900, 
1300 and 2200 counts in the </, r and i bands, respectively. Using (1) it is obvious 
that the readout noise is negligible, and signals at the sky level are background 
noise limited. Typical counts around the center of galaxies are 10000-20000 ADU, 
making the readout noise even less significant there. The sky noise that is actually 
measured off the galaxy in the images is slightly higher than suggested by (1) 
due to imperfections in flat flelding and the presence of faint objects and cosmic ray 
events in the region where the noise is measured. 

Keeping the exposure time short for all the observations made it possible to 
obtain more than 150 images (including flats and standard star observations) in 
one night. The centers of some of our brightest galaxies, however, are saturated. 
Saturation was limited to only few pixels around the photometric center of these 
galaxies. We introduced a scheme to repair these regions and to calibrate these 
galaxies (see the next section). 

2.2 Lowell 

Eighty two galaxies were observed with the 1.1 meter telescope of the Lowell 
Observatory on several nights between 1989 March 24 and April 4. A camera with a 
thinned, back-illuminated 320 X 512 pixel RCA CCD was in the Cassegrain focus. 



signal in electrons in a given pixel): 




(1) 
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A 2:1 focal reducer gave an f/4 focal ratio, yielding a scale of 1"35 pixel""^ . The field 
of view is about 7' X 11'. For more details, see Guhathakurta and Tyson (1989). 

We used the three-band R and / photometric system, which was 

introduced to take advantage of silicon CCDs (GuUixson et al., 1994). Only Bj 
and R images were obtained of each galaxy. These filters are centered at 450 and 
650 nm. Frei and Gunn (1994) give system response curves representing the system 
as we used it. 

Typical exposure times were between 100 and 600 s, depending on the galaxy 
core brightness. This yielded sky levels between 100 and 1000 counts (the sky in 
the R band being typically twice as high as in the B j band for identical exposure 
times). Short text exposures were first obtained of each galaxy in each filter in 
order to avoid core saturation. The readout noise of the CCD was 90 electron, and 
the gain of the instrument was 11 electron/ADU. From (1) we conclude that in 
the background the readout noise and the sky noise are about equally important. 
Typical counts in central regions of the galaxies reach above 10000 ADU, so there 
the shot noise dominates over the readout noise. 

Both observing runs took place in spring. Since we tried to avoid galaxies close 
to the Galactic equator, all the galaxies are located in the Northern Galactic cap. 
The positions of the galaxies in the catalog are plotted in Figure 1 in equatorial 
coordinates. We used a Hammer-Aitoff projection which produces an equal area 
map of the entire sky. The Galactic plane is shown with a thick curve for orientation. 
The CCD images are oriented such that North is up, and East is to the left. 



3. Data Processing 



Debiasing and flat flelding the images of galaxies observed at Palomar was 
straightforward. We used the data in the extended horizontal register of the CCD 
to obtain the bias level, which was subtracted from the image. Twilight sky flats 
were obtained in all three passbands during the night of observations. These were 
used to create a median flat which was used to flat fleld the images. 

Besides subtracting the bias level, we had to use a more complicated flattening 
procedure for the images observed at Lowell. Because of the glass substrate of 
the RCA CCD, the night sky emission lines caused Fabry-Perot fringing in the R 
band (Bj is not affected since there are no strong atmospheric emission lines in this 
band). Because of fringing, we treated the R and the Bj images differently. Another 
problem was that observations took place in both bright and dark conditions, since 
the observing run was longer than a week. The color of the background sky is 
close to the color of the twilight flat in bright conditions, but can be quite different 
if the Moon is down. Consequently, we handled images taken in bright and dark 
conditions differently. Some of these techniques were used by Tyson (1986). 

For the Bj images, if the background was bright, we divided by a twilight 
flat. Because the colors are similar, this was sufficient to obtain good results. If 
the sky was dark, after dividing by the twilight flat we also subtracted a carefully 
scaled twilight flat to account for the color difference. For the R frames, where 
the fringing amplitude can be as large as 15% of the night sky level, in addition to 
the steps performed for the Bj frames, we also subtracted a scaled "fringe frame". 
Twilight flats have a much lower fringe amplitude than even the bright sky images. 
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We scaled the fringe frame to remove fringes as best we could from the twilight flat, 
then used it to flat-fleld the object exposures. We constructed the fringe frame by 
imaging relatively sparse parts of the sky for 300 s at least flve times each night, 
and by averaging these frames by median flltering to minimize contributions from 
individual stars and galaxies. 

3.1 Core Saturation in Palomar Images 
A few pixels around the photometric centers of some of the galaxies observed at 
Palomar were saturated. Compared to the total area of the flelds, only a very small 
portion of the images were affected, accounting for only a negligible part of the 
total light of the galaxy. We modeled the light distribution around the photometric 
centers - L(x^y) - with a seven-parameter, two dimensional Gaussian distribution 
of the form: 

L{x,y) = A + Bexp[-Ci{x-af-C2iy-bf-Cs{x-a){y-b)] , (2) 

where a and b are the coordinates of the photometric center of the given galaxy, 
X and y are the rectangular coordinates in the image, A, 5, Ci, C2, and C3 are 
constants. These 5 constants and the coordinates of the center of the galaxy are to 
be fltted. 

First, to preserve integer representation of the images, we divided all pixel 
values by 2 in the entire image if the center was saturated; in no case did this 
compromise the sampling of the noise histogram. Next, we used the information 
in the image just around the saturated center to flt these parameters, and then 
replaced all saturated pixels with analytical values from (2). Since the few images 
which are saturated were overexposed only slightly, dividing by 2 was enough in all 
cases to ensure that the fltted centers are lower than the maximum value, 32767, 
that can be represented by short signed integers in the pixels. We made a clear 



indication whether the images were saturated (and fixed by this method) in the file 
headers (see Section 5). 

The fitting algorithm was the multidimensional downhill simplex method of 
Nelder and Mead (1965). The actual code is adopted from Press et a/.'s (1988) 
amoeba routine. Crucial for the success of the fit is the selection of sensible starting 
values of the parameters to be fitted. The estimated center of the galaxy (a and 
b) was obtained as a user input; the background sky level (A) and the cross-term 
(C3) were set to 0. The widths in x and y dimensions (Ci and C2) were estimated 
from the full width at half maximum (FWHM) of the point spread function (PSF) 
representing the image. The peak value (B) was set twice as high as the saturation 
level (32767). With these starting values the iteration converged quickly and gave 
good results. 

Eleven of the 31 spiral galaxies observed at Palomar were saturated at the 
photometric center, eight of them in all bands, three of them only in two bands. 
All saturated images were fitted with the above method and their saturated pixels 
replaced. Different extents of saturation in the different passbands resulted in the 
Gaussian distribution being fitted to different regions of the same galaxy. Because 
of this, colors may be unreliable in the very central regions of these galaxies. 

3.2 Star Removal 

We used an empirical two dimensional PSF to fit and remove the foreground 
stars from the images. The program we developed for this purpose first identifies 
stars off the galaxy to be used for constructing the PSF. Second, it finds those 
objects which are likely to be foreground stars (versus HII regions, bright stars 
belonging to the galaxy, etc) and removes them using a PSF fit. Third, the residuals 
left in the image are repaired cosmetically. Although the program does not perform 



all necessary steps automatically, it was very useful for the purpose of cleaning these 
images. Details of this procedure are to be found in a separate paper (Frei, 1995b). 

On several occasions the area occupied by the galaxy was much smaller than 
the total area of the original image, and we chose to extract square regions out of 
the original images that contained the galaxy, which saved a considerable amount of 
storage space. The CCD chip used at Lowell was not square, making it necessary to 
trim to obtain square images. We decided to construct the catalog using only square 
images of a few standard sizes. Offsets (the coordinates of the lower left corner of 
the trimmed image in the original image) and the coordinates of the photometric 
center of the galaxies in the final, clipped images were recorded in the file headers 
(see Section 5). 

We also calibrated all the galaxies photometrically, but not without difficulty. 
Conditions during several nights of observations at Lowell, and part of the observing 
night at Palomar were non-photometric. Moreover, observations of standard stars 
were not available for all nights at Lowell. In a separate paper (Frei and Gunn, 
1994), we derived transformation relations for colors of galaxies among five different 
photometric systems. These were used to calculate magnitudes in the given filter 
systems from data available in the Third Reference Catalog of Bright Galaxies (de 
Vaucouleurs et al., 1991, RC3). For those data from Palomar, both the calculated 
magnitudes and the magnitudes obtained from standard-star-calibrations are given 
in the file headers (there were seven standard stars observed during the night 
of observations). The Lowell data contain only magnitudes calculated via color 
transformations of the RC3 data. 

We summed up total counts corresponding to the galaxy in a circle centered on 
the photometric center of the galaxy (the x and y coordinates of which is recorded 
in the file header). The radius of this circle was half of the horizontal (or vertical. 
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since all images are square) size of the image (which is also recorded in the file 
header). This aperture is large enough to effectively make our magnitude the total 
magnitude of the galaxy. We subtracted the contribution of the sky using the total 
area of the circle and the sky level recorded in the file header. If the photometric 
center of the galaxy does not coincide with the center of the image, we used only the 
part of the circle which is inside the image for the sum, and calculated the area of 
this part for correct sky subtraction. With all the necessary information available 
from the FITS header keys, these steps can be reproduced with ease. Starting from 
this total count representing the galaxy, and the calibrated magnitude in the given 
band (derived via color transformation of RC3 data), we obtained the count which 
would represent magnitude in a pixel; this is stored in the file header. 



4. Properties of Galaxies in the Catalog 



All of the galaxies in our catalog are listed in RC3. We recorded some of the 
properties found in RC3 for these galaxies in the file headers (such as morphological 
type, total B magnitude, B — V color, etc). Here we present distributions of some 
of the most important properties, to statistically describe the catalog. In the next 
section, we will give detailed information about all the parameters included in the 
file headers. Although we developed techniques — some of them new — to derive 
similar properties and to automatically process the images for classification, the 
aim of this paper is to report the data, not our reduction technique. Thus we 
chose to present parameters from RC3, a well known source. We will present the 
properties we obtained with our methods for these galaxies elsewhere (Frei, 1995a, 
in preparation). 

The revised morphological type on the "Handbuch der Physik" system is used 
in RC3 (de Vaucouleurs, 1959, 1963; de Vaucouleurs eA al., 1991). The histogram 
of the numerical type T of the galaxies in our catalog is in Fig. 2. Among the 
113 galaxies, one is of type 10 (non-Magellanic irregular). The numerical type (T) 
is 90 for this galaxy in RC3. We omitted this galaxy from the histogram. The 
distribution of morphological types is well suited to our classification task. We have 
many more spirals than ellipticals. Deriving parameters of spirals (based on the fine 
structure in the disk) is more complicated than finding representative parameters 
of ellipticals (where fine structure, although well studied, is not usually considered 
for classification purposes). 

Total Bt magnitudes are derived from photoelectric aperture photometry and 
from surface photometry in RC3. The histogram in Fig. 3 shows the distribution 



of this parameter. Although we cahbrated our images in their respective bands, 
there is not a single band in which observations are available for all the galaxies (we 
used different photometric systems at Palomar and at Lowell). Consequently, the 
easiest way to compare brightnesses is to use data from RC3. As a rough measure 
of the distances of the galaxies, we plotted the histogram of the radial velocities 
corrected to the reference frame of the cosmic microwave background radiation (u3k ) 
in Fig. 4. Isophotal major axis diameters (-D25) listed in RC3 were obtained from 
photographic and CCD surface photometry and from photoelectric growth curves. 
We produced a histogram of these sizes in Fig. 5, where the diameters are expressed 
in arcseconds. The majority of our galaxies are 4-6 arcmin in diameter, some of 
them even larger; our catalog contains some of the brightest, largest, most nearby 
galaxies around the Northern Galactic cap. 

RC3 lists isophotal axis ratios for the galaxies, both ellipticals and spirals. We 
calculated inclination angles (I) from these data according to 



where R25 is the decimal logarithm of the axis ratio (as it is actually given in 
RC3, obtained at an isophotal surface brightness level jib = 25 mag arcsec"^). As 
usual, 0° represents face-on galaxies, while 90° corresponds to edge-on orientation. 
We plotted a histogram of the inclination angles for our spiral and lenticular 
galaxies (99 in total) in Fig. 6, since the axis isophotal ratio is not directly related 
to an inclination angle for ellipticals. The apparent lack of face-on galaxies is 
understandable since spiral galaxies, even if face-on, will have non-circular isophotes 
at any brightness level due to the uneven surface brightness distribution of light in 
the disk (spiral arms). Edge-on, or almost edge-on galaxies will have isophotes 
which appear "fatter" than the orientation of the disk would indicate, due to finite 




(3) 



disk thickness and due to the presence of the bulge. This wiU resuh in the measured 
inchnation angle (from (3)) being somewhat (systematically) larger than the "true" 
disk inclination. 

In Fig. 7 we plot the distribution of position angles. The standard convention 
for measuring position angles is from the North eastward between 0° and 180°, as 
it is given in RC3. It is reassuring that they are evenly distributed (there is no 
unreasonable selection effect in our catalog). The apparent excess of galaxies in 
the first bin (0 position angle) is because the position angle was artificially set to 
for any nearly face-on spirals in RC3 (for those with inclination angles less than 
% 40°). 

Some of the most important properties of the galaxies are given in Table 1. 
The first column is the NGC number of the catalog galaxy. The observatory 
where the images were obtained is identified in column 2. All Palomar galaxies 
are in three bands (three separate images are available, in </, r and z), while all 
Lowell galaxies are in two bands (two images, in Bj and R). The numerical 
revised morphological type T from RC3 is in column 3; these data are used to 
construct the histogram in Fig. 2. Columns 4 and 5 are right ascensions and 
declinations as recorded in the file headers at the time of observations. Columns 
6 to 10 are the total Bt magnitudes, radial velocities corrected to the reference 
frame of the microwave background radiation (i^sk), major axis isophotal diameters 
(-D25), inclination angles, and position angles, respectively. These data are all from 
RC3, and they were used for Figs. 3-7. 
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5. FITS header keys 



We recorded useful information in the header keys in the FITS files of the 
images. There are several keys among those we used which are standard FITS keys 
(for the description of the FITS format and the standard keys see Wells et al., 1981; 
or "Definition of the Flexible Image Transport System (FITS)", from NOST). The 
standard keys that we populated with data are: SIMPLE, BITPIX, MXIS, MXISl, 
MXIS2, OBJECT, TELESCOP, EPOCH, BSCALE, BZERO, BUNIT, DATAMAX, DATAMIN and 
DATE-OBS. 

There are several keys that we added to describe the data as best as we could. 
All necessary information is in Table 2. The first column contains the name of the 
header key, the second column has the type of the variable, and the last column 
gives a short description. The first six keys (EXPOSURE to TIME) were recorded 
during the observations. We obtained SKY and SKYSIG from the data at the time of 
calibration (the sky level is the median of all pixels in the image). The photometric 
zeropoints DNATO_ST and DNATO_BV are calculated from observations of standard 
stars (not available for the Lowell data) and by using data from RC3, respectively. 
We recorded the offsets of the trimmed, smaller image in the original, full size 
image and the photometric center of the galaxy in the trimmed, square image. The 
parameter SATURATE is used to tell whether the galaxy center was saturated and 
repaired. The last seven header keys are populated with data from RC3. These are 
well described in RC3 and in the previous section. We gave the location (column 
and line) of the actual data entry in the RC3 table. 
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6. Conclusion 



During 1989-1990, a substantial set of CCD surface photometry of nearby, 
well resolved galaxies were collected from Palomar and Lowell. The data were 
used to construct a catalog of images in FITS format, containing a uniform set 
of information in the file headers. The catalog is suitable for our ongoing galaxy 
classification project (Frei, 1995a, in preparation), and we believe it can also be used 
for other purposes. We collected and presented information about the observations 
and various data processing steps in the previous sections. The description of FITS 
header keys and brief statistics of key properties of the catalog galaxies are also 
included to aid the prospective user of these data. 

To our knowledge, no data similar to that presented here are available to the 
general public today. Although we plan to increase the number of galaxies in 
the catalog, we think that the current set of images is already useful for various 
applications, and thus it should be made accessible to anyone interested. The 
image files are on our "anonymous" FTP site at "astro.princeton.edu". The data 
are in "frei/Galaxies" . Additional information is in the README file in the same 
directory. The data may also be accessed via World Wide Web (WWW) at URL 
"ftp: / / astro.princeton.edu/frei/Galaxies" . 



- 20 - 



7. Acknowledgements 



It is our pleasure to acknowledge the great number of people who participated 
in the observations and data reductions in the past 5 years. Their work made it 
possible to construct the catalog as presented here (we obviously try to avoid saying 
"to complete", since such a large set of data can be improved and extended almost 
indefinitely). Students at the Department of Astrophysical Sciences, Princeton 
University, did observations at Lowell as part of a student project. We thank Jill 
Knapp for organizing the project. We are indebted to several people at Lowell: 
Jay Gallagher, who made it possible to use the observatory for the project, Craig 
GuUixson, who helped with the telescope, Neal Hartsough (University of Arizona) 
and Rick Wenk (AT&T Bell Labs) who built the software system. We are greatly 
indebted to the former students who took part in the data acquisition and reduction: 
Arif Babul, Renyue Cen, Charles Gammie, Bin Gao, Lynne Hillenbrand, Neil 
Katz, Man Hoi Lee, Joanna Lees, Kevin Long, Shude Mao, Jordi Miralda-Escude, 
Changbom Park, Duncan Walsh and David Weinberg. We acknowledge very useful 
discussions with Robert Lupton and Neil Tyson, who also helped to revise the 
manuscript. The patience of the Palomar time allocation committee, who continued 
to give time for this project through five completely clouded-out years, is also much 
appreciated. 

This research was supported in part by NSF through grant no. AST-9100121. 
Z. F. acknowledges the support of OTKA through grant no. F4491. 



- 21 - 



REFERENCES 

Boroson, T. 1981, ApJS, 46, 177 
Burstein, D. 1979a, ApJS, 41, 435 
Burstein, D. 1979b, ApJ, 234, 435 

CoUess, M., Burstein, D., Wegner, G, Saglia, R. P., McMahan, R., Davies, R. L., 

Bertschinger, E., & Baggley, G. 1993, MNRAS, 262, 475 
Davis, L. E., Cawson, M., Davies, R. L., & lUingworth, G. 1985, AJ, 90, 169 
de Vaucouleurs, G. 1959, Handbuch der Physik, 53, 275 
de Vaucouleurs, G. 1963, ApJS, 8, 31 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Buta, R., Paturel, G., & 
Fouque, P. 1991, Third Reference Catalog of Bright Galaxies (Springer- 
Verlag, New York) 

Freeman, K. C. 1970, ApJ, 160, 811 

Frei, Z. 1995a, ApJ, to be submitted 

Frei, Z. 1995b, PASP, submitted 

Frei, Z., & Gunn, J. E. 1994, AJ, 108, 1476 

Guhathakurta, P., & Tyson, J. A. 1989, ApJ, 346, 773 

GuUixson, C. A., Boeshaar, P. C., Tyson, J. A., & Seitzer, P. 1994, ApJS (in press) 
Gunn, J. E., & Knapp, G. R. 1993, in Sky Surveys: Protostars to Protogalaxies, 

ASP Conference Series 43, edited by Soifer, B. T., p. 267 
Gunn, J. E., & Westphal, J. A. 1981, Proc. SPIE, 290, 16 
Han, M. 1992, ApJS, 81, 35 

Hester, J. J., & Kulkarni, S. R. 1989, ApJ, 340, 362 
Kent, S. M. 1984, ApJ, 56, 105 



- 22 - 

Kent, S. M. 1985, ApJ, 59, 115 
King, I. R. 1962, AJ, 67, 471 
King, I. R. 1965, AJ, 70, 376 
King, I. R. 1966, AJ, 71, 64 
King, I. R. 1978, ApJ, 222, 1 

Kodaira, K., Okamura, S., & Ichikawa, S. 1990, Photometric Atlas of Northern 

Bright Galaxies (University of Tokyo Press, Tokyo) 
Kormendy, J. 1977a, ApJ, 214, 359 
Kormendy, J. 1977b, ApJ, 217, 406 
Kormendy, J. 1977c, ApJ, 218, 333 
Lauer, T. R. 1985, ApJS, 57, 473 

Nelder, J. A., & Mead, R. 1965, Computer Journal, 7, 308 

NOST: NASA/Science Office of Standards and Technology 1993, Definition of the 
Flexible Image Transport System (FITS), Code 633.2 (NASA Goddard 
Space Flight Center, Greenbelt, MD) 

Okamura, S., Kodaira, K., & Watanabe, M. 1984, ApJ, 280, 7 

Press, H. W., Flannery, B. P., Teukolsky, S. A., & Vetterhng, W. T. 1988, Numerical 
Recipes in C (Cambridge University Press, Cambridge) 

Sandage, A. 1961, The Hubble Atlas of Galaxies (Carnegie Institution of 
Washington, Washington, D. C.) 

Sandage, A., & Bedke, J. 1988, Atlas of Galaxies (NASA, Washington, D. C.) 

Sandage, A., & Tammann, G. A. 1981, A Revised Shapley-Ames Catalog of Bright 
Galaxies (Carnegie Institution of Washington, Washington, D. C.) 

Schneider, D. P., Gunn, J. E., & Hoessel, J. G. 1983, ApJ, 264, 337 

Simien, F., & de Vaucouleurs, G. 1986, ApJ, 302, 564 

Thuan, T. X., & Gunn, J. E. 1976, PASP, 88, 543 



- 23 - 

TuUy, R. B., & Fisher, J. R. 1987, Nearby Galaxies Atlas (Cambridge University 

Press, Cambridge) 
Tyson, J. A. 1986, J. Opt. Soc Am. A, 3, 2131 

Wade, R. A., Hoessel, J. G., Elias, J. H., & Huchra, J. P. 1979, PASP, 91, 35 

Watanabe, M. 1983, Ann. Tokyo Astr. Obs., 2nd ser., 19, 121 

Watanabe, M., Kodaira, K., & Okamura, S. 1982, ApJS, 50, 1 

Watanabe, M., Kodaira, K., & Okamura, S. 1985, ApJ, 292, 72 

Wells, D. C, Greisen, E. W., & Harten, R. H. 1981, A&AS, 44, 363 

Whitmore, B. C, & Kirshner, R. P. 1981, ApJ, 250, 43 

Yoshizawa, M., & Wakamatsu, K. 1975, A&A, 44, 363 



- 24 - 



FIGURE CAPTIONS 

Figure 1. Positions of the 113 galaxies in the catalog on the sky. An equal-area 
Hammer-AitofF projection was used. The center is at 12h RA and 0° declination. 
The Galactic plane is shown with the thick curve. All galaxies are around the 
Northern Galactic cap, and avoid the Galactic plane. 

Figure 2. Distribution of the numerical morphological types T of the galaxies 
in the catalog. T is obtained from RC3 for all the galaxies. Types —5 and —4 are 
ellipticals, —3 to —1 are lenticulars, to 9 are spirals and the rest are irregulars. 

Figure 3. Distribution of the total Bt magnitudes of the galaxies in the 
catalog. Bt is from RC3 for all the galaxies. 

Figure 4. Distribution of the velocities, corrected to the reference frame of 
the cosmic microwave background radiation, of the galaxies in the catalog, vsk is 
from RC3 for all the galaxies. 

Figure 5. Major axis isophotal diameters of the galaxies in the catalog. These 
diameters, at a surface brightness level jib = 25 mag arcsec"^, are found in RC3 
for all the galaxies. 

Figure 6. Distribution of inclination angles of the spiral, lenticular and 
irregular galaxies in the catalog. The data are found in RC3 for all the galaxies. 
The apparent lack of face-on and edge-on galaxies are due to the techniques used 
to determine the inclination angles. See Section 4 of the text for an explanation. 

Figure 7. Distribution of position angles of all galaxies in the catalog. All 
data are from RC3. The high number of galaxies in the first bin is due to the fact 
that all nearly face-on spiral galaxies were assigned position angle in RC3. 
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